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We describe direct imaging measurements of a colloidal sphere's diffusion between two parallel 
surfaces. The dynamics of this deceptively simple hydrodynamically coupled system have proved 
difficult to analyze. Comparison with approximate formulations of a confined sphere's hydrodynamic 
mobility reveals good agreement with both a leading-order superposition approximation as well as 
a more general all-images stokeslet analysis. 



Stationary surfaces modify the flow field set up by 
a moving particle, thereby increasing the particle's hy- 
drodynamic drag. Calculating this influence remains a 
vexing problem in all but the simplest geometries be- 
cause complicated boundary conditions usually render 
the problem intractiblc. For example, Faxen derived a 
single sphere's hydrodynamic coupling to a rigid planar 
surface as early as 1924. Adding a second parallel 
wall, however, is so much more challenging that a defini- 
tive formulation is not yet available. This is particularly 
unsatisfactory in that many biologically and industrially 
relevant processes are governed by particles' dynamics in 
confined geometries. 

This Letter describes measurements of a single col- 
loidal sphere's diffusion through water in a slit pore 
formed by two parallel glass surfaces. We use optical 
tweezers Q to position the sphere reproducibly within 
the sample volume, a slow steady Poisseuille flow to es- 
tablish its position between the walls, and high-resolution 
digital video microscopy || to track its motions in the 
plane. By positioning a test sphere at selected heights 
within a slit pore, releasing it and tracking its motions, 
we are able to measure its height-dependent hydrody- 
namic coupling to parallel bounding surfaces. These mea- 
surements agree quantitatively with predictions based on 
the stokeslet approximation, a tool which is particularly 
useful for describing many-body hydrodynamic interac- 
tions in colloidal suspensions . 

Previous imaging Q and light-scattering Q studies 
probed confined spheres' dynamics averaged over the 
slit pore's thickness and only indirectly addressed how 
a sphere's mobility changes as it moves relative to con- 
fining walls. A very recent study combined digital video 
microscopy with optical tweezer manipulation to mea- 
sure a sphere's mobility near the midplane of a slit pore 
as a function of the slit pore's width ||. Its results cast 
serious doubt on a recently proposed theory for confined 
Brownian motion 0, but left open questions regarding 
a sphere's dynamics at other, less symmetric configura- 
tions. 

The present measurements were performed on a sin- 



gle polystyrene sulfate microsphere from a suspension 
of spheres a = 1.006 ± 0.010 fim. in radius (Catalog 
No. 4204A, Duke Scientific) dispersed in an aqueous so- 
lution of 2 mM NaCl at T = 300.66 ± 0.07°K. The sam- 
ple was cleaned before resuspension by extensive dialysis 
against deionized water and then infiltrated into a slit 
pore of area 2 cm x 1 cm and thickness H — 8 /im cre- 
ated by sealing the edges of a #1 coverslip to a glass 
microscope slide with uv cured adhesive (Norland Type 
88). Access to the sample volume was provided by two 
glass tubes bonded to holes drilled through the slide at 
either end of the longest dimension. All glass surfaces 
were cleaned thoroughly before assembly to ensure uni- 
form surface properties IJ. Although both the sphere 
and the glass walls develop large surface charges when 
immersed in water, the suspension's high ionic strength 
reduced the Debye-Huckel screening length to 7 nm and 
thus minimized electrostatic interactions. Surface sep- 
arations were always large enough that van der Waals 
attractions were negligible [ |To| |. 

The sample was mounted on the stage of an Olympus 
IMT-2 inverted optical microscope and imaged with a 
lOOx NA 1.4 oil immersion objective lens. Images cap- 
tured with an attached CCD camera were recorded on 
a JVC BRU-S822DXU computer-controlled SVHS video 
deck before being digitized with a MuTech MV-1350 
frame grabber. Computerized analysis || of the result- 
ing sequence of digitized images yielded measurements of 
the colloidal sphere's position in the microscope's focal 
plane with 20 nm spatial resolution at 1 /60 sec intervals 

§• 

A colloidal sphere diffuses through a Newtonian fluid 
according to 

(Arf(r))=2A(f)r, (1) 

where Ar^(r) is its displacement in the ?-th dimension 
over time r and Di(r) is the associated diffusion coeffi- 
cient. Angle brackets indicate an ensemble average. The 
fluctuation-dissipation theorem relates Di(r) to the drag 
force, 7; (r) Vi , the sphere experiences as it moves past 
point r with speed v^. Di{f) — fc^T '/ji(r). The drag 
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FIG. 1. (a) Measured probability distribution for displace- 
ments along the direction of flow at r = 1/60 sec and 4/60 sec 
for 2 = 4 /jm. Solid curves are fits to Eq. ^ for the width 
and displacement of the distribution, (b) Drift of the distri- 
butions' centers along and perpendicular to the direction of 
the imposed flow, (c) Evolution of the mean-square widths of 
P(n,r\z), fit to Eq. (Eh for D^h). 



coefficient, 7^, emerges from the solution to the Stokes 
equation describing flow at low Reynolds numbers, sub- 
ject to the surfaces' no-slip boundary conditions. 

For example, a sphere of radius a suspended in an 
unbounded fluid of viscosity r\ has an isotropic and 
translation-invariant drag coefficient, 70 = Girrja. The 
corresponding free self-diffusion coefficient for the spheres 
in the present study is D = 0.259 ± 0.04 /xm 2 /sec. A 
sphere passing between two rigid parallel walls, on the 
other hand, experiences a drag which depends both on its 
position z in the slit and also on its direction of motion. 
Measuring this dependence requires the ability either to 
track a sphere's motion in three dimensions, or else to 
position the sphere reproducibly within the slit. 

We adopted the second approach, using an optical 
tweezer to place a single sphere in the microscope's fo- 
cal plane, releasing it for t = 5/60 sec to measure its 
diffusivity, and then retrapping it. The optical tweezer 
was formed by directing a collimated beam of laser light 
(100 mW at 532 nm) through the objective lens' back 
aperture. Optical gradient forces exerted by the tightly 
converging beam localized the sphere near the focal point 
despite radiation pressure and random thermal forces. 
Suddenly deflecting the beam onto a beam block with a 
galvanometer-driven mirror extinguished the optical trap 
and freed the particle to diffuse ||. Residual heating of 
roughly 1°C due to steady-state optical absorption re- 
laxed in a few microseconds through thermal diffusion, 
and did not affect the sphere's dynamics on the time scale 
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FIG. 2. The drift rate's dependence on initial height in the 
slit pore. Solid curves result from least squares fits to Eq. (Q). 
Shaded regions are forbidden by sphere- wall contact, given the 
fit estimates for the wall locations. 



of our data collection H,[ll) . The galvanometer drive was 
synchronized to the video camera's sync signal to ensure 
that data acquisition began at a reproducible interval af- 
ter the optical trap was extinguished. 

We adjusted the focal plane's height in 1.0 /mi in- 
crements to sample the diffusivity's dependence on the 
sphere's height in the slit pore. The resulting uncertainty 
in z due to the sphere's out-of-plane diffusion was smaller 
than \J2D§t = 0.2 /im, and so was small compared with 
both the sphere's diameter and the separation between 
walls. 

We determined the absolute positions of the glass walls 
by observing the influence of a steady Poisseuille flow on 
the particle's motion. The advected sphere attained a 
height-dependent drift velocity 
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where zq is the position of the lower wall. This drift 
caused the probability distribution for finding the sphere 
at position r to shift as it broadened with time: 
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where (^(r)) = Vi(z)r is the mean displacement in the 
j-th direction after time r given the steady flow in that 
direction at height z. 

At each tweezer height z, we collected 3,000 se- 
quences of 5 video fields, yielding 12,000 measurements 
of P(x,t\z) and P(?/, t\z) over four 1/60 second intervals. 
Typical examples at the midplane of the slit pore appear 
in Fig. |(a). The peaks' drifts, plotted in Fig. g(b), yield 
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FIG. 3. Height dependence of the in-plane diffusion coef- 
ficient for a Brownian sphere in a slit pore. The solid curve 
results from the stokeslet approximation and is indistinguish- 
able from the linear superposition approximation's prediction. 
Dashed curves show the range of predicted values due to the 
uncertainty in the sphere's radius. Squares show Faxen's pre- 
dictions from Eq. (o). 



components of the drift velocity, Vi(z), which appear in 
Fig. ^. Fitting Vi(z) to Eq. (||) for ifyj, H, and Zo, yields 
H = 8.16 ± 0.22 /im and fixes the lower wall's absolute 
position to within ±0.13 /xm. Errors primarily reflect the 
experimental uncertainty of ±0.25 /im in the increment 
between tweezer locations. Without loss of generality, we 
set zq = and measure starting heights h = z — zq from 
the lower glass surface. 

The evolution of the distributions' widths over time, 
shown in Fig. 0(c), can be interpreted with Eq. (|l|) to 
obtain the drag-corrected self-diffusion coefficients Di(h) 
for each height. Results are summarized in Fig. ||. As 
expected, the in-plane diffusion coefficients measured for 
motions along and transverse to the imposed flow agree 
to within 5% at all heights. The sphere's diffusivity is 
strongly suppressed near either wall, and falls well below 
Do even along the midplane. 

Faxen's result for the drag on a sphere near a single 
wall provides a good starting point for interpreting these 
observations. His result for the in-plane drag, 71 (/i), 

71(h) 16 h \h 3 J { ' 

has been verified to great accuracy using photonic force 
microscopy in a thick sample cell [H. Oseen suggested 
in 1927 that the drag due to two parallel walls could be 
approximated by a sum of two such single-wall contribu- 
tions 0: 

12(h) w 70 + [11(h) - 70] + [li(H -h)- 70]. 

(5) 



Even though this linear superposition approximation vi- 
olates boundary conditions at both bounding surfaces, 
the corresponding prediction for the in-plane diffusivity, 
D 2 (h) = kgT/j 2 (h), agrees quite well with our measure- 
ments, with no adjustable parameters, as can be seen in 
Fig. [| Higher-order corrections ffl to 71 (h) yield negligi- 
bly small corrections to Z^ft.) and so are not presented. 

Faucheux and Libchaber adopted the linear superposi- 
tion approximation, averaged the resulting in-plane dif- 
fusivity over h, and obtained reasonably good agreement 
with /i-averaged measurements for various sphere radii 
B. Lin, Yu and Rice further demonstrated linear super- 
position's accuracy at h = H/2 by measuring a sphere's 
in-plane and out-of-plane diffusivity near the midplanes 
of slit pores of various widths || . 

Regardless of its empirical success, the linear superpo- 
sition approximation fails to satisfy boundary conditions 
and so cannot yield accurate predictions under all cir- 
cumstances. Lobry and Ostrowsky attempted to remedy 
its shortcomings by accounting for flows' reflections off 
bounding surfaces |Q . Considering such corrections up to 
third order, they obtained results which compare rather 
poorly with measurements by Lin, Yu and Rice Q . 

Faxen managed to calculate two-wall drag coefficients 
for the particularly symmetric arrangements h = H/A 
and h = H/2 Q. The in-plane results, 

{1 - 0.6526^+0.1475^ 
-0.131^-0.0644^, ft=f 
1 - 1.004 1 + 0.418 £ (6) 
+0.21^-0.169^, ft=f 

appear as squares in Fig. |^ and agree well with both 
Oseen's approximation and also with our measurements. 
Unfortunately, Faxen's method does not work for other 
configurations. 

Another approach, due to Blake [Q, takes advantage 
of an analogy between hydrodynamics and electrostatics. 
Rather than attempting to satisfy no-slip boundary con- 
ditions at a bounding wall directly, Blake introduced the 
notion of a hydrodynamic image whose flow field exactly 
cancels the particle's on the boundary. Uniqueness of the 
Stokes equation's solutions guarantees that the image so- 
lution also solves the original problem. 

Blake's treatment yields the Green's function for a par- 
ticle's flow field bounded by a single surface. Account- 
ing for a second parallel surface involves adding not only 
the particle's hydrodynamic image in the new surface, 
but also an infinite series of images of the images. Liron 
and Mochon [fl5f summed this series explicitly for a point 
force acting on the confined fluid, derived the additional 
terms needed to cancel residual flows on the surfaces, and 
thereby obtained the Green's function, G(hz — r, H), for 
the flow at f due to a disturbance a height h above the 
lower wall in a slit pore of width H. Integrating their 
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result over the sphere's surface would yield an accurate 
though unwieldy expression for j(h). 

A much simpler solution can be obtained when far field 
contributions dominate the flows at the bounding sur- 
faces. This is the case when the sphere is much smaller 
than the gap between the walls, a -C H . In this limit, 
the flow field due to the sphere's hydrodynamic images 
is well approximated by the image flow for a point force. 
Thus, the drag on a sphere centered a height h above the 
lower wall may be approximated by 



1 - 1 (h)^ lQ 1 + [G(0,H)-G s (0)}, 



where Gs(0) is the singular contribution due to a unit 
point force at the sphere's location, so that the term 
in square brackets is the flow at the sphere's location 
due to the images alone. More generally, Gs(r) is the 
Green's function in an unbounded fluid and is known as 
a stokeslet. Eq. (^), then, is an example of a stokeslet 
approximation. Evaluating Eq. (^) using the expressions 
for G(r) and G s (r) in Equations (5), (15), (32) and (33) 
of Ref. pa] yields the in-plane diffusion coefficient 
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is the contribution due to a stokeslet at the sphere's po- 
sition together with an infinite series of image stokeslets, 
and 
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Xrj 2 sinh 2 [A(l - rj)\ + coth(A) sinh 2 (A?7) - 



rf) sinh(Ary) — r\ sinh[A(2 — i])] 
sinh(A) 



(10) 



enforces the no-flow boundary conditions at the walls. 
The numerically evaluated result appears as the solid 
curve in Fig. || 

Remarkably, the predictions of Eqs. (8), (JsJ) and 
(10) are indistinguishable from Oseen's superposition ap- 
proximation for our experimental conditions, and both 
agree quantitatively with Faxen's fifth-order results, 
Eq. (^) . Further comparison with Eq. (^) reveals that the 
stokeslet approximation becomes increasingly accurate 
for larger wall separations, while the linear superposition 
approximation fares less well. For smaller separations, on 



the other hand, the stokeslet approximation becomes less 
accurate. Linear superposition performs surprisingly well 
at small separations, by contrast, particularly if higher 
order corrections to Eq. (^) are included. 

Despite the apparent complexity of Eqs. (8), (0) and 
(10) when compared with the linear superposition ap- 
proximation, stokeslet analysis has the appeal of scala- 
bility. Under conditions for which the stokeslet approxi- 
mation is valid, contributions to the diffusivity tensor for 
a system of spheres can be combined in a straightforward 
manner 0], as in Eq. (^). We recently demonstrated this 
approach's utility in the comparatively simple case of two 
spheres near one wall for which linear superposition yields 
inaccurate predictions (jsj. Further demonstrating its ac- 
curacy for the much more challenging case of two-wall 
confinement lends confidence in its potential for system- 
atically treating many-body colloidal hydrodynamics. 
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